First and second sound in cylindrically trapped gases 
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We investigate the propagation of density and temperature waves in a cylindrically trapped gas 
with radial harmonic confinement. Starting from two-fluid hydrodynamic theory we derive effective 
ID equations for the chemical potential and the temperature which explicitly account for the effects 
of viscosity and thermal conductivity. Differently from quantum fluids confined by rigid walls, 
the harmonic confinement allows for the propagation of both first and second sound in the long 
wave length limit. We provide quantitative predictions for the two sound velocities of a superfluid 
Fermi gas at unitarity. For shorter wave-lengths we discover a new surprising class of excitations 
continuously spread over a finite interval of frequencies. This results in a non-dissipative damping 
in the response function which is analytically calculated in the limiting case of a classical ideal gas. 
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Superfluids are known to exhibit, in addition to usual 
sound, an additional mode, called second sound, where 
the superfluid and normal components move with oppo- 
site phase PP. In a weakly compressible fluid, like su- 
perfluid helium, second sound reduces to a temperature 
wave, leaving the total density practically unaffected. 
In helium its velocity was systematically measured as a 
function of temperature, providing a high precision deter- 
mination of the superfluid density. First measurements of 
second sound were recently reported also in dilute Bosc- 
Einstein condensed gases confined in elongated traps [2] . 
In this Letter we show that the propagation of sound in a 
nonuniform trapped gas exhibits new interesting features 
both in the superfluid and in the normal phase. We will 
focus on highly elongated configurations where the radial 
harmonic trapping provides the relevant non-uniformity. 
These configurations are well suited for the experimental 
excitation and detection of sound waves [3H5]- We will 
consider the usual hydrodynamic regime I <C A where I is 
the mean free path and A is the wavelength of the sound 
wave. We will also consider the condition of strong radial 
confinement R± <C A where R± is the radial size of the 
sample. 

The above conditions are compatible with two differ- 
ent regimes, depending on the ratio between the viscous 
penetration depth 8 = yjrj/ p n ui and the radius R± where 
r\ is the shear viscosity coefficient, p n is the normal den- 
sity and uj is the frequency of the sound wave. Let us 
first consider the case of a uniform fluid confined by the 
hard walls of a tube. If 8 <C viscosity plays an unim- 
portant role and sound propagates similarly to the case 
of bulk matter. If instead 8 3> R± , the viscosity imposes 
the uniformity of the normal velocity field as a function 
of the radial coordinate. Since friction further requires 
that the normal velocity be zero on the walls, the nor- 
mal part of the fluid cannot move at all along the tube 
and only the superfluid can. This mode is known as 4-th 
sound [6 . In the case of a harmonically trapped gas the 



two regimes exhibit new features whose investigation is 
the object of the present paper. 

Our analysis is based on the two-fluid hydrodynamic 
equations (HD) [7] 

a 4 p+V-j = 0,m9 t v s = -V(/i + C/), (1) 
dth + d,P + pdiU/m = d k {r,T lk ), (2) 
d t s + V • (sv n ) = V • (kVT/T) , (3) 

which reduce, above Tc, to the classical equations of hy- 
drodynamics. In the above equations j = p s v s + p n v„ 
is the current density, p s and p n are the superfluid and 
normal mass densities for a fluid with total mass density 
p = p s + p n , v s and v„ are the corresponding veloc- 
ity fields, T lk = (d k v ni + d l v nk - 28ikdjV n j/3), s is the 
entropy density, while P and p are, respectively, the lo- 
cal pressure and chemical potential. U is the external 
trapping potential which will be assumed of axial har- 
monic form: U — m{u>\r\ + uj 2 z 2 )/2 with r\ = x 2 + y 2 
and uj z <C Since we are interested only in the lin- 
ear solutions, in the above equations we have omitted 
terms quadratic in the velocity. Finally r\ and k are, re- 
spectively, the shear viscosity and the thermal conductiv- 
ity (we have not considered bulk viscosities terms which 
give smaller contributions). In a uniform fluid the ef- 
fect of viscosity and thermal conductivity is irrelevant 
in the long wave length limit, providing only higher or- 
der corrections to the dispersion law. In the presence of 
confinement their effect is instead crucial. This opens 
new perspectives to explore the behavior of viscosity and 
thermal conductivity in interacting trapped gases. The 
two-fluid hydrodynamic equations, in the dissipationless 
regime (77 — n — 0), have been already the object of nu- 
merical investigations [8HT0] in the presence of harmonic 
confinement. These works, so far limited to isotropic 
trapping, have provided predictions for the discretized 
frequencies in the superfluid phase with special emphasis 
to the hybridization effects between the first and second 
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sound solutions. 

In this Letter we are interested in the solutions of the 
hydrodynamic equations satisfying the condition ui <C 
cjj_. As a consequence of the tight radial confinement, 
Eq. (J2J for the current implies the important condition 
X7±P + pVj_U = of mechanical equilibrium along the 
radial direction. Violation of this condition would in fact 
result in frequencies of the order of ujj_ . The tight radial 
confinement also implies that the radial component of the 
velocity field must be much smaller than the longitudinal 
one. 

Low frequency regime. In the presence of harmonic 
trapping the condition S 3> R± of large viscous penetra- 
tion depth is equivalent to requiring the low frequency 
condition lo <C lo\t where r is a typical collisional time. 

Excluding v s from Eqs. |2]) and ([T]) and using the ther- 
modynamic identity dP — sdT + pdp/m one can write 
the equation for the relevant z component of the velocity 
field of the normal component in the form 

mp n d t v z n + p n d z (Sfi) + msd z (ST) = mV • (r?V<) , (4) 

where we have ignored terms containing the small radial 
components of the velocity field. 

The presence of viscosity in Eq. Q results in the inde- 
pendence of v z n on the radial coordinate r± . In fact vio- 
lation of such a behavior would be incompatible with the 
low frequency condition us <C io\r. It is worth noticing 
that, differently from the case of the tube with hard walls 
discussed in the introduction, for harmonic trapping the 
uniformity of the velocity field does not stop the motion 
of the normal component. Analogously, the presence of 
thermal conductivity in Eq. Q for the entropy results in 
the independence of the temperature fluctuations on rj_ 
in the same low frequency limit This in turns im- 

plies that also the fluctuations of the chemical potential 
will be independent of the radial coordinate. This follows 
from the radial mechanical equilibrium condition and the 
use of the thermodynamic identity dP — sdT + pdp/m. 
Thus in the low frequency limit both the fluctuations ST 
and Sp are independent of the radial coordinates. 

It is now convenient to integrate radially both the 
equations of continuity for the density and for the en- 
tropy ([l]|3| as well as Eq. Q for v z nl profiting of the fact 
that the large terms containing the derivatives with re- 
spect to r_L vanish due to Gauss theorem, and that those 
containing the second derivative with respect to z can be 
ignored since they give rise to higher order corrections 
in the dispersion law. Integration of Eq.Q then gives 
—m{p n )dtv z l = m(s)d z 5T + (p n )d z 5p, where the symbol 
(...) = J ...dxdy stands for the integral on the transverse 
variable rj_. The above equation for v* and the analog 
equation mdtv z g = —d z 5p for the superfluid component 
can be used to derive the following coupled equations for 



the averaged density and entropy fluctuations: 

- d 2 t 5{ P ) + d z ((s)d z ST) + d z (( P )d z Sp) jm = 0, (5) 
-dfS{s) + d z ((s) 2 /(p n )d z ST) + d z ((s)d z Sp) /m = 0.(6) 

One can further write S(p) = d^(p)Sp + dx(p)ST and 
S(s) = d ll (s)Sp + 8t(s)ST The spatial dependence of the 
thermodynamic quantities entering the above equations 
is fixed by the trapping potential U through the local 
density approximation equilibrium condition p(p, T) + 
U = po, with pq fixed by the normalization condition for 
p. In the presence of axial trapping (lo z 7^ 0), Eqs. ( 5|6 ) 
allow for the calculation of the low frequency discretized 
states of the system. In the following we will assume 
lj z = (cylindrical geometry) in order to calculate the 
velocity of the sound waves propagating along the axial 
direction. The solutions for Sp and ST have the t and z 
dependence cx e l (9^- wt ) with lo — cq, yielding the most 
relevant equation 

c 4 [md^p)d T (s) - (d T (p)) 2 } + c 2 [2d T (p)( S )- 
( P )d T (s) - md,(p)(s) 2 /(p n )} + (p s )( S ) 2 /(Pn) = 0, (7) 

which generalizes the well known Landau equation [TJ [7] 
for the first and second sound velocities to the case of a 
cylindrically trapped gas. 

Above Tc, where the superfluid density vanishes, 
Eq.|7]) admits only one solution with velocity different 
from zero. In this case the value of c can be more con- 
veniently calculated by direct radial integration of the 
HD Eqs. ([lj3]). In the low frequency limit the equations 
of continuity for the density and the entropy yield the 
ID isoentropic condition dt{(s) / (p)) = and one imme- 
diately finds the result c 2 = {d{P) /d{p))( s )/( p ) , which 
generalizes the most famous expression for the adiabatic 
sound velocity of uniform matter as well as the zero tem- 
perature result mc 2 — (d fl (p) / (p)) -1 derived in [T5]. At 
unitarity, as well as for the ideal gas, where sj p depends 
on the combination pT -3 / 2 , the new adiabaticity con- 
dition takes the form (p)T~ 5 / 2 = const. A non trivial 
implication concerns the sound velocity in the ideal clas- 
sical limit where one finds c 2 = (7/5) (T/m) rather than 
the value (5/3) (T/m) holding in uniform gases. Here and 
in the following we put the Boltzmann constant ks = 1- 

Before discussing the solutions of Eq. ^ in the pres- 
ence of harmonic trapping, it is useful to calculate the 
dynamic response function of the system to an external 
potential of the form Ae t ( qz ~ ut K The response function 
is determined by the density fluctuations induced by the 
external potential according to S(p) = \x(q,uj)e lt - qz ~ u ' t K 
Its calculation is important because it provides informa- 
tion on the actual possibility of exciting the two sound 
waves using a density probe. The inclusion of the exter- 
nal perturbation affects both the equation for the current 
and the equation for v s . Calculations yield the following 
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result for the imaginary part of the response function 



Imx(q,cj) = -</>)-{WiCj[5(a;-ciq) + £(u; + ciq)] 



(27rT 5 / 2 /w^) J^°^ T v(u)du and analogously for the other 
quantities. In conclusion all the coefficients of the disper- 
sion relation (JtJ) can be calculated as a function of T and 
of the dimensionless variable u$ = po/T. In practice, in 
order to determine the function H(u) we have used a fit 
to the experimental data of [5D] above the critical tem- 
perature Tc and assumed that, below Tc, the pressure is 
constant as a function of T, following the analysis in the 
same reference. The results for the sound velocities ci 
and c 2 and for the relative weights W2/W1 [see Q] are 
reported in Fig. [T]for the unitary Fermi g clS 0jS el function 
of T/T c in the relevant interval 0.5T C < T < T c , where 
Tc ~ 0.197V, having defined the Fermi temperature for 

our cylindrical geometry as T F = (^(|^) 3/ Vl(p)) 2/5 
[22] . In Fig. [2] we show the ratio between the relative 
density (S(p)/(p)) and temperature (ST/T) fluctuations 
in the two modes. The ratio provides a physical insight 
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+ W 2 u}[S(u-c 2 q)+5(LJ + c 2 q)]}, (8) c 



where the weights Wi and W 2 obey the relationships 
mcfWx + mc 2 2 W 2 = 1 and Wi + W 2 = (d^p) / (p)) T 
determined, respectively, by the /- and the isothermal 
compressibility sum rules. The knowledge of Im% is rele- 
vant for experiments based on the propagation of density 
pulses and can be measured directly in Bragg scattering 
experiments p~3l I14j . 

We are now ready to calculate the values for the two 
sound velocities and the relative weights in the density re- 
sponse. We will specialize to the case of the unitary Fermi 
gas [TS] since in this case the collisional regime, needed to 
apply hydrodynamic theory, is more easily achieved due 
to the large value of the scattering length. At unitarity 
the effects of the interactions in the thermodynamic func- 
tions can be expressed in a universal way |16j in terms of 
the dimensionless parameter u = p/T. Thus the pressure 
can be written as P(p,T) = T 5 ^ 2 H(u) where the func- 
tion H(u) can be determined through microscopic many- 
body calculations [17H19] or, in some ranges of temper- 
ature, directly extracted from experiment [20) . All the 
thermodynamic functions can be expressed in terms of 
the function H . For example the density takes the form 
p(p,T) = mT z / 2 v(u) with v(u) = H'(u), the entropy 
density is given by s = 5P/2T — up/m etc.. Only the 
superfluid mass density p s — mT 3 ' 2 v s (u) requires the 
knowledge of another independent function for which we 
use the regularized form introduced in [21) . 

The radially integrated quantities entering the disper- 
sion relation ([7| can also be easily calculated. Since 
in the presence of radial trapping the chemical poten- 
tial has the radial dependence p(r) = po — mw|r}/2 
the integration over dxdy can be usefully transformed 
into integrals over u. For example, we find (p)(po,T) — 
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FIG. 1: First (green) and second (red) sound calculated us- 
ing Eq.([7|, as a function of T/Tc, normalized to the Fermi 
velocity vf = \/2Tf /m. The inset reports the ratio W2/W1 
between the weights of the two sounds modes in Imx [see Eq. 



on the nature of the modes. In the first sound solution 
(high velocity mode) the relative variations of density are 
larger than the ones of temperature and have the same 
sign. The opposite happens for second sound. The value 
of the ratio W 2 /W\ is smaller than the value found for a 
uniform fluid [13] and this is due to the reduction of the 
superfluid component as one moves from the symmetry 
axis of the trap. We have also checked that our predic- 
tions are not very sensitive to the fitting procedure used 
to calculate the thermodynamic functions. 

High frequency regime. We have so far discussed the 
solutions of the hydrodynamic equations in the low fre- 
quency limit. In the opposite limit u 3> us\t of large 
frequencies viscosity and thermal conductivity can be ig- 
nored :23J. In this limit, still compatible with the condi- 
tion wCwi, the solutions are also affected by the pres- 
ence of the radial confinement, but in a different way. In 
particular the fluctuations of the temperature and of the 
chemical potential are no longer independent of the ra- 
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FIG. 2: Ratio between the relative temperature (ST/T) and 
density {S{p)/{p}) fluctuations in the two sounds as a function 
of T/T c (see Fig. [l]for the notation). 
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dial coordinate, except in the superfluid region [24 . This 
limit is still compatible with the usual hydrodynamical 
condition lot <C 1 only if uj±t <C 1, a condition which is 
rather severe from the experimental point of view. 

An explicit and instructive solution in the dissipa- 
tionless regime is available above Tq where the super- 
fluid component is absent. We will consider the extreme 
regime of the ideal gas where all the thermodynamic 
quantities are known explicitly. In this limit the classical 
equations of hydrodynamics can be recast in the form of 
an equation for the velocity field [25 : 



function in the region of the continuum 



muj v = — TVfV • vl 

3 1 i 



V[v- VZ7] 



-[V-v]VC/. (9) 



Solutions of this equations for 3D harmonic traps were 
considered in [26] and [27] ■ The advantage of using ID 
cylindrical configurations is due to the availability of 
plane wave solutions along the z direction, providing new 
physical insight on the effects of trapping on sound propa- 
gation. It was found in [35] that an exact solution of this 
equation in the presence of cylindrical radial trapping 
(U = mui]_r\/2) is given by v z cx / 5 e lqz , v± — 0, where 



j^_r^_/T, and is characterized by the adiabatic 



velocity of sound mc 2 = 5T/3. This solution is quite dif- 
ferent from the one holding in the low frequency regime 
discussed in the first part of the paper where v z does not 
depend on the radial coordinate. Also the temperature 
and the chemical potential fluctuations exhibit a strong 
radial dependence as a consequence of the radial trap- 
ping and the absence of thermal conductivity and vis- 
cosity. From the explicit knowledge of the velocity field 
and the orthogonality condition J pvf^ • v^dxdy = if 
LJi ^ ujj for the solutions of the HD equations [9] , it is pos- 
sible to calculate the contribution of the adiabatic sound 
to the dynamic response function. Using the expression 
p oc e - ^ / 2 for the density distribution at equilibrium one 
obtains, after straightforward algebra, the result 



Hp) 



9 m q 2 c 2 



(10) 



showing that the adiabatic mode does not exhaust the 
/-sum rule x( w — ^ °°) = ~9 2 (p)/ TOtj2 ; differently 
from the case of a uniform gas, and hence revealing 
that additional solutions of the hydrodynamic Eq. ([9| 
should exist. This equation actually admits an addi- 
tional class of low frequency solutions lying in the con- 
tinuum and consequently exhibiting damping. On the 
u — q plane these solutions occupy the region w/q < 
^/24/25c a = Co- The velocity field is given by the 
expressions v z = Ae^l 4 [crcos(a-£ 2 ) — sin(cr£ 2 )/20] e lqz , 
\/ c^q 2 /lo 2 — 1/4 is a continuous param- 



where a 



Be? I 4 sin^ 2 



7f, with B 



eter, and vj_ - 

—iAy/3/20 [{c 2 a q 2 - ijj 2 )/qc a uj±\ . We present here the 
result for the imaginary part of the dynamic response 



ImXc(?,w) 



4 (p) ujy/clq 2 - uj 2 



3 toCq 



c 2 <? 2 



(11) 



It is easy to show that inclusion of the continuum allows 
for a complete fulfillment of both the compressibility and 
the /-sum rules. The function Imx c exhibits a sharp 
maximum at ui/q 0.98co- Thus the new excitation can 
be considered as the analog of second sound above Tc- It 
is worth emphasizing that the damping associated with 
the continuum is not a consequence of dissipative effects 
as happens in uniform fluids above Tc, but a result of 
the energy leakage from the gas due to the increase of 
the velocity field on large distances. The existence of 
a continuous spectrum of excitations and the value of 
the limiting velocity cq are universal properties holding 
also at lower temperatures, being fixed by the asymptotic 
behavior at large r± where the gas is classical and ideal. 
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